Nonequilibrium effective vector potential due to pseudospin exchange in graphene 
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We show that exchange interactions in two-dimensional electron gases out of equilibrium can gen- 
erate a fictitious vector potential with intriguing signatures in interference and Hall measurements. 
Detailed predictions are made for graphene, where the effect is enhanced by pseudospin exchange. 
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By driving a system out of equilibrium one typically 
corrupts its quantum coherence. As a result many quan- 
tum effects observed in electronic structures, such as the 
Aharonov-Bohm interference and the Kondo effect are 
suppressed as one applies a voltage bias [J, 13, S] ■ In rare 
instances, however, deviations from thermal equilibrium 
modify quantum effects in more interesting ways, or give 
rise to new coherent dynamics, such as the AC- Josephson 
effect In this Letter we predict a novel effect in the 
latter category: a fictitious vector potential induced by 
nonequilibrium conditions. The effect is particularly ro- 
bust and intuitive in graphene [^, [li 0] ' where the con- 
duction electrons obey an equation that resembles the 
relativistic Dirac equation [8|. In that equation an or- 
bital degree of freedom that is usually referred to as the 
'pseudospin' plays the role of the conventional electron 
spin in the Dirac equation. This new quantum number 
has profound consequences, such as 'Klein tunneling', al- 
lowing electrons to pass arbitrarily high potential barriers 
, and a velocity renormalization through exchange in- 
teractions [lOj . Here we show that in the nonequilibrium 
state created by a voltage bias the same pseudospin ex- 
change interactions induce an effective vector potential. 
Unlike the fictitious vector potentials created by lattice 
defects or distortions in graphene If, 13, l3| this effect 
docs not cancel between the two bandstructure 'valleys', 
but it has directly observable consequences. We discuss 
two of those signatures: i) interference currents oscillat- 
ing as a function of the current density in the material 
and ii) a 'Hall' voltage that reveals the fictitious magnetic 
field implied by the predicted vector potential in curved 
conductors. 

The origin of this effect is best understood within the 
Dirac model for electrons in graphene at low energies. 
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where — (cr^jjay) is a vector of Pauli matrices in 
pseudospin space, A is the magnetic vector potential, 
7 = ± the valley index, p the electron momentum, and 
V, — e, c are the Fermi velocity, the electron charge, and 
the speed of light, respectively. Eq. (P) imphes that in 
each valley electrons have a definite projection of their 
pseudospin onto their velocity. When a bias voltage is ap- 
plied that drives an electrical current, such that electrons 



have a preferred direction of motion, a graphene sheet 
therefore carries an excess pseudospin. Consequently one 
expects a pseudospin-dependent exchange energy of elec- 
trons, described by an effective Hamiltonian with a vector 
potential ^4'^'' in addition to A in Eq. ([T]). 

Exchange vector potential: To quantify the effect that 
is anticipated by the above simple Hartree-Fock argu- 
ment we first consider a clean, infinite sheet of extrinsic 
graphene, with Fermi energy ep > 0, at A = 0. The 
effects of electron-electron interactions on single-particle 
quantities, as considered here, are captured by the elec- 
tron self-energy E. We evaluate S in the so-called GqW- 
approximation, resulting in a 'screened self-energy'. The 
'Fock'-contribution to the retarded self-energy then is 
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Here, the noninteracting retarded and Keldysh [f 4:] Green 
functions G*^"-'^ and G^"^^, respectively, have an implicit 
matrix structure in pseudospin space and they are diag- 
onal in the spin and valley indices <t and 7. The retarded 
and the Keldysh components of the RPA interaction po- 
tential = vc/{l - wcn^) and = HKjVF'^p^ re- 
spectively, are found from the unscreened Coulomb in- 
teraction wc(q) = 27rrsD/|q| at interaction parameter Ts 
and the polarizability H. The bare Keldysh Green func- 
tion G^*^-*^ is characterized by occupation numbers fsk of 
single-electron states ijja-ysk- Here, k are the electron mo- 
menta and s = f , ~f for the electrons in the conduction 
band and the valence band, respectively. Through G*^°^^ 
also and depend on fsk- We consider a state 
with current flow in the direction we — (cos 0, sin 0), 
fsk = B[eF -I- e(sfe • ■W0)eV - e^fc], with Q{x) = 1 for 
a; > and 0(a;) = for x < 0. Here, Sgk are the single- 
particle energies, V is the voltage that drives the cur- 
rent, and we take the limit of temperature T = 0, requir- 
ing kT <ti eV in practice. We assume weak nonequilib- 
rium, eV <C ep , such that we may expand in small devia- 
tions 6 fsk — V{dfsk/dV\v=Q) of the occupation numbers 
from their equilibrium values. We find the leading order 
(in eV/ep) nonequilibrium contribution to the self-energy 
S'E,^^ ~ V{dY^^' /dV\v=o) at small Tg by expanding the 



(T7 



2 



G(")K of Eq. © in Sfsk, 
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X M^i^(e-e,fc,p-fc). (3) 

Here, k — kw^ and takes its equilibrium form, 
obtained in Refs. [3, At eV <ti ep, 6fsk is non- 

vanishing only if \ssk ~ £f| ^ ^F- For all quanti- 
ties computed below S is needed, in the same limit, 
at \e — EfI <C ep and thus \e — e^fcl ^ £f in Eq. ([3]). 
Since the typical momentum transfer p — A; in Eq. ([3]) 
is of order fcp, we may substitute W^{e — esk,P — fc) w 
W^ {0,p-k) = 2TTrsv/{\p-k\+qsc) into Eq. ©, with the 
screening wavevector qsc = 4rsfcF at Fermi wavevector 
fcp [H, Ell • The pseudospin-diagonal components of the 
real part Re ^E^'^"^ shift the chemical potential. The off- 
diagonal components of the hermitian part of (5E^(°) have 
a p-wave contribution that implies the velocity renormal- 
ization analyzed in Ref. fi^ and an s-wave part that acts 
as an effective vector potential A^^, 



= -— Rc / — Tr<T^JS^f (fcFii>^,eF) 
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(we set ?i = 1), which, for the above fsk, evaluates to 
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At rg <C 1 we have C ^ Ar^ ln(7r/8rs)/7r and in the range 
of 0.15 < Ts < 2.4 typical for graphene [l^ C is of 
order unity (C ~ 0.38 at ^ 0.15 and C ~ 0.89 at 
Ts = 2.4). Note that the valley-dependence of the effect 
indeed drops out, as advertised [l8|. 

Expanding the and the of Eq. ^ in 5 fsk re- 
sults in corrections to C. These corrections take the form 
of integrals over higher powers of than the leading 
contribution and they are thus negligible at rg ^ 1, when 
the GoW-approximation is rehable We therefore 

employ our 'screened Hartree-Fock' approximation Eq. 
([3]) in the remainder of this Letter. The Keldysh self- 
energy E^ and the anti-hermitian part of E^ account for 
inelastic relaxation. These inelastic processes are negligi- 
ble in the geometries we study below: graphene ribbons 
with a finite length L that are ballistic, that is they have 
long inelastic and elastic mean free paths Zin and Id, re- 
spectively, Zin, lp\ ^ L, and they are well-coupled to elec- 
tron reservoirs [20|. Neglecting relaxation by phonons 
(T = 0), we find from Eq. ^ for a transport electron 



with momentum k = kw^f, in the above nonequilibrium 
state fsk at <C 1 and eV <C rsSp, 



6e'^{3eV -2de) 



37re| 



5i Vl^"'^ + W2|+4rg 



(6) 




'A\l+eV 



FIG. 1; Two ballistic graphene ribbons attached to large 
leads (black bars) , biased by a voltage V such that a current 
flows through the upper ribbon. Electrons tunnel between the 
ribbons at Xi and X2 with matrix element A, driven by a volt- 
age 5V. The amplitudes for electron propagation between xi 
and X2 along the upper ribbon with exchange vector poten- 
tial A**^ interfere with those along the lower ribbon, causing 
oscillations of the tunnel current as a function of V. 



with Se ^ Esk - ep, valid at \Q - <i)± 7r/2| > ey/ep 20]. 

An effective vector potential A!^^ is generally expected 
in electron systems where the exchange energy between 
two electrons depends on their relative motion. This 



requires gsc 



< kT, 



which is typically fulfilled in n-type 

In 
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GaAs, but not in p-type GaAs or Si-MOSFETs 
graphene the effect persists into the strong screening limit 
qsc ^ k-p (at Ts ^ 1) through pseudospin exchange. 

Interference: We now turn to a discussion of signatures 
of the predicted effective vector potential. We first con- 
sider two tunnel-coupled, ballistic graphene ribbons, as 
shown in_Fig.[Tl where A^^ can be observed interferomet- 
rically ,22,] . The two ribbons are well-coupled to electron 
reservoirs at their ends and they touch at x\ and X2, 
which creates a weak tunnel coupling A between them. 
For simplicity we assume A to be identical at x\ and xi 
and constant over the width W of the ribbons. We take 
the semiclassical limit fcpVF ^ 1 and allow for different 
path lengths Ax" and Ax' between the points of tunnel- 
ing along the upper and the lower ribbon, respectively. A 
voltage y ^ ep/e that leaves the mean electric potential 
invariant (A/i -ey/2 + ^6% - {eV/2Y - ep in Fig. [J) 
drives a current / through the upper ribbon. According 
to Eq. ^ an A'^^ results that shifts the phase of an elec- 
tron traveling from xi to X2 along the upper ribbon by 
$ = CeyAx"/87rw. The tunneling current /tun between 
the ribbons, driven by a small voltage 5V , is expected to 
oscillate as a function of We accordingly find 
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for the differential tunnel conductance Gtun — 
dItun/d5V\sv=o at generic ep and |Aa;" — Ax^\ <C Ax" 
pot . Here ~ l/kpW depends on the Fermi wavevec- 
tor of the lower ribbon kp. We conclude that in the ge- 
ometry of Fig. [1] the predicted exchange vector potential 
A"^^ has the very distinctive signature that the tunnel 
conductance Gtun between the two ribbons oscillates as a 
function of the current flowing through the upper ribbon. 

The visibility of the predicted oscillations is quanti- 
fied by the interference contrast C — Gint/Gdir. We ex- 
tract the 'direct' conductance Gdir and the oscillatory 
'interference' conductance Gint from Eq. ([7]) as the first 
two and the third terms in curly brackets, respectively. 
We find G ~ 1 for_fcFAa;^ > %/]V (/x e {u,l}), but 
N. Larger path length differences 



Ax 



:$> VN suppress the interfer- 
on 



/cfIAx" - Ax' I < 
fcpAa;^ > A:f|Ax" 

ence contrast, G ~ VTV/fepI Ax" — Ax'jlnA'^. Above we 
assumed tunneling matrix elements A that are constant 
over the width of the ribbons, so that electrons tunnel 
only between identical transverse modes. To estimate the 
influence of inhomogeneities of A we compute the tunnel- 
ing conductance also for random tunnel matrix elements 
with identical ensemble averages for scattering between 
any two of the transverse modes. We find G ~ and 
we expect that G suffers similar or weaker suppressions 
for most realistic forms of tunnel contacts. 

Violations of our limit Ax'^ ^ kn suppress G expo- 
nentially. According to Eq. ([6]) the condition Ax^ ^ hn 
is met at low voltages. Demanding at the same time 
that V be large enough to produce one entire oscilla- 
tion of Gtun imposes a lower bound on the path lengths 
fcpAx'' > 327r2rs/^C"^''^ [11 ■ Note that Eq. ^ assumes 
a two-dimensional electron gas, W 3> v/eV. In the one- 
dimensional limit, realizable in carbon nanotubes, quasi- 
particles with well-defined energies cease to exist and the 
above predictions do not apply. Inelastic or elastic scat- 
tering that induces a potential drop AV along the rib- 
bons may also cause an electric Aharonov-Bohm effect. 
In order for the resulting oscillations not to obscure the 
predicted effect it has to be assured that AV^ ^ v/eAx^. 

Hall effect: In a curved conductor, cut out of a sheet of 
graphene, the current direction wg is not constant any- 
more. We thus expect an exchange vector potential A^^ 
similar to Eq. but space-dependent. The effective 
magnetic field that would generically result should be 
observable by the Hall effect that it causes. To confirm 
this scenario we consider a segment of a graphene an- 
nulus with inner radius R, radial width W, and angular 
width (3 (see Fig^). We again address the semiclassical 
limit kpW ^ 1 [231 and omit the details of the calcula- 
tion that will be presented elsewhere [2^ . The segment is 
well-coupled to electron reservoirs at its angular bound- 
aries, such that electrical current flows in the angular 
direction and we have periodic boundary conditions at 
the angular edges. The wavefunctions with these bound- 
ary conditions obtain from the solutions of Eq. ([T]) in 




FIG. 2: Graphene annulus of width W and inner radius 
R, contacted by electron reservoirs (solid black bars) whose 
electrochemical potentials differ by eV. A current / is flowing 
that produces a transverse voltage Kr by a Hall effect in the 
(fictitious) exchange magnetic field B'^''. 



the full annulus, where angular momentum I is a good 
quantum number. Additionally the single-particle states 
i^a-ysni are characterized by a radial quantum number n 
and they have eigenenergies Ssni- At Tg = a voltage bias 
between the two reservoirs in Fig.[2]creates an occupation 
of states fsni = Q[eF+e{sl)eV-es„i] (T = 0). Substitut- 
ing the corresponding G*^°^^ into Eq. in our 'screened 
Hartree-Fock' approximation corresponding to Eq. ([2]), 
and extracting A"^^ through Eq. (|4]) yields an exchange 
vector potential A^^{rWip) — Cic/v){V/8n)w^_^_^/2 as 
suggested by Eq. ([5|), pointing in angular direction. 
The above A'^^ implies an effective magnetic field 
= jA'^'^l/r perpendicular to the annulus. The cor- 
responding Lorentz force modifies the charge density p, 
found from the Green functions implied by Eq. ([3]), 
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Here, and follow from Eq. ([T]) with A = A^ 

and in the semiclassical limit we obtain 
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The equilibrium contribution po is r-independent 
The r-dependent contribution to p, Eq. Q, is screened 
by a charge rearrangement through a transverse voltage 
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+ 0{e-'Vye'p), (10) 



restoring charge neutrality. Eq. (jlOp predicts a transverse 
'Hair voltage generated by an electrical current in the 
absence of a real magnetic field. It is a second peculiar 
signature of the exchange vector potential A®^. 

Eq. (fTO|) crucially relies on our assumption that the 
annulus is ballistic, Im ^ L, where L = max{/3i?, W} 
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and = niin{Zci, Zin}-In the opposite limit of a diffusive 
annulus, L ^ l^, the current density decays as 1/r, im- 
plying A^^ cx it'(^_|_7r/2 Z*" with B'^^ — and thus V^^.^ = 0. 
With moderately strong backscattering, ~ L, A^^ still 
produces a transverse voltage. Backscattering, however, 
induces a voltage drop along the annulus that generi- 
cally produces an additional transverse voltage 1/^'^'^^'*. 
The linear in V contribution to V^^^^^'t can be elimi- 
nated by measuring Vn- — [Vtr{V) + Vtri—V)]. Any 
energy-dependence of the scattering, however, will add 
a quadratic in V contribution to V^^^'^^*' that competes 
with V^^^. Typically backscattering i n g raphene has 
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an energy-dependence on the scale ep [25j, |26 
It follows that V'tf'^" is negligible, V^f''" < ^^t^, if 
L/ln, < ln(l -I- W/R). Experimentally V^f ^" can be 
distinguished from V^^.^ by varying ep at eV <C £f- If 
the measured Vtr is due to energy-dependent scatter- 
ing, a variation of ep by c?/i ~ eV causes a relative 

change of Ftr by AV;«^catt/y^scatt ^ (y^scattj-y^l 



V'scatt 
tr 



tr V ;i£F+ti/i ~ 

1. In contrast, the relative 
dfi/ep ^ 1 if Vtr is due to 



change is only AV^f/V^,^ 
the exchange interaction. 

A second effect competing with V^™ is the Hall effect 
in the (real) magnetic field produced by the current flow- 
ing in the annulus. In the limit W/R ^ 1 the relative 
magnitude of the resulting transverse voltage t/j.™^^" is 
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, 8i?tan(/3/4), 3 
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up to corrections of order {W/ R)\n{W/ R). For the pa- 
rameters W = 0.2R, (3 = 7r/6, and = 0.15 we find for 
instance V;^^8"/V;°'^ « 4 • 10"'' kpW . For typical sample 
dimensions of w 1/im this effect is thus negligible up 
to the largest fcp ~ 1 nm. Also the smallness of VJ.™''^" 
can be verified experimentally by its distinctive depen- 
dence on £y. In this case V^™''^" is independent of ep, 



^^magn 



^^magn ^ compared to AV^^^/V^^ ~ d/.t/£p. 
Both competing effects discussed above can be ruled out 
by simultaneously varying £p and V . 

Conclusions: We have predicted a novel nonequilib- 
rium effect in coherent two-dimensional electronic sys- 
tems. The effect is due to electron-electron interactions 
and it generates a fictitious vector potential that has 
striking signatures: interference currents that oscillate 
as a function of the current density in the material and a 
'Hall' voltage in zero magnetic field. We have made de- 
tailed predictions for graphene, where the effect is partic- 
ularly robust owing to the pseudospin degree of freedom 
of the conduction electrons. 

The author thanks C. W. J. Beenakker, P. W. Brouwer, 
R. Mani, Y. V. Nazarov, and A. Zangwill for discussions. 
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